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ABSTRACT 

The stellar population, metallicity distribution and ionized gas in the elliptical galaxies 
NGC 6868 and NGC 5903 are investigated in this paper by means of long-slit spec- 
troscopy and stellar population synthesis. Lick indices in both galaxies present a neg- 
ative gradient indicating an overabundance of Fe, Mg, Na and TiO in the central 
parts with respect to the external regions. We found that M§2 correlates both with 
FelA527o an d FeiA5335, suggesting that these elements probably underwent the same 
enrichment process in NGC 6868. However, only a marginal correlation of Mgi and 
Fe±A527o occurs in NGC 5903. The lack of correlation between computed galaxy mass 
and the M<?2 gradient suggests that these elliptical galaxies were formed by merger 
events. The stellar population synthesis shows the presence of at least two populations 
with ages of 13 Gyr and 5 Gyr old in both galaxies. 

We have estimated the metallicity of the galaxies using SSP models. The central 
region of NGC 6868 (\R\ < 0.5 kpc) presents a deficiency of alpha elements with respect 
to iron and solar metallicity. The external parts, present a roughly uniform distribution 
of [a/Fe] ratios and mctallicitics ranging from [Z/Zq] = —0.33 and solar. A similar 
conclusion applies to NGC 5903. 

Concerning the emitting gas conspicuously detected in NGC 6868, we test three hy- 
potheses as ionizing source: an H II region, post-AGB stars and an Active Galactic Nu- 
cleus (AGN). Diagnostic diagrams involving the ratios [NII]\ 65S 4/ Ha, [OI] A6300 /Ha 
and [SII]x67i7,3i/Ha, indicate that values measured in the central region of NGC 6868 
are typical of LINERs. Together with the stellar population synthesis, this result sug- 
gests that the main source of gas ionization in NGC 6868 is non-thermal, produced by a 
low-luminosity AGN, probably with some contribution of shocks to explain ionization 
at distances of ~ 3.5 kpc from the nucleus. 

Key words: Elliptical galaxies: NGC 6868 and NGC 5903 - fundamentals parameters 



1 INTRODUCTION 

A question still open to debate is how elliptical galaxies 
are formed. Three scenarios have been proposed to ad- 
dress this issue: (i) monolithic co llapse of gas clouds linke d 
to high energy diss ipation rates dArimotp fe Yoshii 1 1 19871 ): 
(ii) galaxy merger IjBurkert fe Naabll2000T ): and (Hi) hybrid 
models that involve an early monolithic collapse fo llowed by 
merger events at later times l|Ogando et aUT 2005b A conse- 
quence of the monolithic collapse (associated to a large-scale 
star formation) is that a conspicuous metallicity gradient is 
expe cted to establish along the full range of gala ctic radius 
fe.g. lCarollo et al.lll993al : IChiosi fe Carraroll2002l ). However, 
if merger events are important, the gradient is expected to 
disappear. A few ellipticals present strong signs of recent 



dynamical perturbation ( Schweizer et al.lll990l ). which may 
provide important clues to their formation processes and 
evolution. 

Important information on the above i ssues has been de - 
rived from metal line-strength indices fe.g. lDavis et al.lll987l ) 
and their radial variation within galaxies. Mg2 line-strength 
distribution in early-type galaxies, for example, can vary 
considerably, ranging from essentially featureless to struc- 
tured profiles showing, e.g., changes of slope possibly as- 
sociated with kinematically decoupled cores, or anomalies 
in the stellar population of some ellipticals (|Carollo et al.l 
ll993bh . 

In addition, the origin of the ionized gas in ellipticals is 
another controversial issue, because until recently it was be- 
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lieved that the interstellar medium (ISM) in these galaxies 
was n ot significant. However, studies by, e.g. Phillips et al.l 
(1986) have shown that more than 50% of the elliptical and 
lenticular gala xies contain signifi c ant am ounts of ionized gas. 
More recently, iMacchetto et all (|l996l ) have traced ionized 
gas in a sample of 73 elliptical and SO galaxies by means 
of narrow-band filters centered on Ha and [Nil]. They de- 
tected emission gas in about 3/4 of the galaxies, with masses 
amounting to 10 3 to 10 5 M Q . 

In the present work we investigate properties of the el- 
liptical galaxies NGC 6868 and NGC 5903, because previous 
studies hav e shown that they contain important int erstel- 
lar medium (|Macchetto et al 1 19961 ; iFerrari et al.l|l999l) have 



complex stellar and gas kinematics (|Caon et al 



2000). Be- 



sides, both ellipticals belong to galaxy groups, which might 
imply past close encounters (and even mergers) with com- 
panions. Potential consequences of such past interactions are 
the triggering of discrete star formation eve nts with stellar 
popu lations differing in metallicity and age l|Worthev et al.l 
1992]) and the presence of large-scale dynamical perturba- 
tions that may affect galaxy morphology by building subsys - 
tems such as a stellar and/or dust disk l|Caon et al.l l2000). 
In this context, NGC 6868 and NGC 5903 can be taken as in- 
teresting candidates where the stellar population and metal- 
licity distribution can be investigated, which are important 
parameters to understand the formation and evolution of 
these galaxies. 

The main goal of the present paper is to investigate, by 
means of long-slit spectroscopy, the origin, star-formation 
history, metallicity and ionized gas distribution in NGC 6868 
and NGC 5903. The analyses will be based mostly on a stel- 
lar population synthesis method (based on templates built 
with star clusters and H II regions) , and the radial distribu- 
tion of absorption and emission lines. 

This paper is organized as follows. In Sect. [2] we describe 
the observations. Lick indices are measuremed are discussed 
in Sect. [3] The stellar population synthesis and results are 
described in Sect. [4] Metallicity and ionized gas are discussed 
in Sect. [5] Discussions are in Sect. [6] 



2 THE TARGET GALAXIES AND 

SPECTROSCOPIC OBSERVATIONS 

NGC 6868 is classified as E3 in RCS^. It is part of the 
GR28 cluster which contains as well the galaxies NGC 6861 , 
NGC6870, NGC6851 and NGC6861D jMaia et al.lll989h . 
Figure [T] shows an R band image of the field containing 
NGC 6868 in two scales (obtained from the NASA/IPAC 
Extragalactic Database - NED). 

The effec t ive rad ius of NGC 6868 is Re = 40" 
( Carollo et al. 1993a). NGC 6868 is a radio source 
(|Savage et al.l 119771 ) and emits in the infrared as indicated 
by 60 and 100 /im IRAS observations. The total luminosity 
in the medium infrared is 15.2 x 10 8 Lq, and the dust mass 
estimated from gra in emission considerations is m 70 Mq 
|Ferrari et alj|2002h . 

The velocity field and velocity dispersion profiles are 




Figure 1. Bottom panel: R image of NGC 6868 showing a scale 
of 12.75' X 12.04'. Top panel: same as before for the region 4.25' X 
4.02' centered on NGC 6868. 



symmetrical about the galaxy center . However, the gas is 
not observed to follow regular orbits l|Zeilinger et al.|[l996l ). 
The latter work shows evidence in several position angles 
that the rotation curve does not remain flat in the outer 
parts but bends back towards a value close to the systemic 
velocity. 

NGC 5903 (morphological type E2 - RC3) forms a 
pair with the elliptical galaxy NGC 5898; both galaxies 
are the brightest members of a small group composed 
of N GC 5903, NGC 5898 and ESO514-G003 l|Maia et all 
1989). NGC 5903 presents a small rotation around the 
semi- major axis dSparks et all 1 19851). Dust is observed in 
a small region f~ 100 pc) (|Ferrari et al]|l999l ). NGC 5903 
presents X-ray emissio n with no evidence of ionized gas 
jO'Sullivan et alj|200ll). The effe ctive radius of NGC 5903 
is Re = 37" l|Carollo etailll993 



1 The Third Reference Catalogue 
Jde Vaucouleurs et al.lll99lh . 



of bright galaxies 



B images (from NED) of 
the field containing NGC 5903 are shown in Fig. [2] 

NGC 6868 and NGC 5903 were observed with the 3.6 m 
ESO telescope (in Chile) equipped with EFOSC1. The ob- 
servational setup included the grism 0150 for the range 
5140 to 6900 A, with a dispersion of 3.4 A pixel -1 , and the 
CCD#2 (512 x 512 Tektronix). This setup was used in or- 
der to maximize the spectral coverage A A 5100 — 6800A and 
obtain the largest possible number of stellar absorption fea- 
tures and gas emission lines wit hout compromisin g the spec- 
tral resolution of 3. 4A pixel -1 (|Caon et alj|2000l ). The spa- 
tial scale of the observational setup was 0.6" pixel -1 with a 
3.1' long slit; the slit width was fixed at 1.5", approximately 
equal to the seeing. 
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Figure 2. Bottom panel: B image of the field of NGC 5903 show- 
ing a scale of 14.16' X 13.37'. Top panel: same as before for the 
region 3.54' X 3.35' centered on NGC 5903. 



At least two spectra were obtained along the position 
angle in order to increase the signal-to-noise ratio and for 
cosmic ray removal. For NGC 6868 three long-slit spectra 
were obtained with 40 min of exposure time; the slit was ori- 
ented along the semi-major axis (PA = 120°). Two 40 min 
spectra were obtained for NGC 5903 with the slit oriented 
along the semi- major axis (PA = 340°). The seeing var- 
ied from 1.2" to 1.8". Several spectrophotometric standard 
stars were observed at the beginning and end of each night 
in order to flux calibrate the spectra. The spectra were 
processed with standard IRAF tasks. The mean bias was 
subtracted from each spectra, which were subsequently di- 
vided by a normalized dome flat-field to remove variations 
in sensitivity. The centroids of the non-saturated lines in 
the He+A comparison spectra were measured in order to 
obtain the wavelength calibration and to map the line cur- 
vature. Galaxy and standard stars spectra were rectified 
and rebinned to a logarithmic scale with a constant step 
of 85kms _1 . The sky contribution for each galaxy spectra 
was determined by taking the median value between two 
windows of 20 pixels wide at both ends of the slit and sub- 
tracted. 

Spectra were extracted in both galaxies for the central 
and 7 regions symmetrically distributed along the north- 
south direction. Pairs of spectra obtained at the same dis- 
tance of the galaxy center were averaged together in or- 
der to increase the signal to noise ratio. We have estimated 
in less than 1% the relative spectrophotometric errors be- 
tween both extractions at each position. The distance of 
each extraction to galaxy center is given in cols. 1 to 3 of 
Table [2] respectively in terms of angular, effective and ab- 
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Figure 3. Rest-frame spectra extracted from NGC 6868, cor- 
rected for foreground reddening. The distance of each extraction 
to the galaxy center is given in arcsec. North and south extrac- 
tions are shown in the left and right panels, respectively. Arbitrary 
constants have been added to the spectra for clarity. Conspicuous 
emission lines (especially H a and [Nil] AA6548, 6584) occur in all 
spectra. 



solute radii. For redshift correction we used the conspicu- 
ous absorption line Nal5S95 as reference. Foreground red- 
de ning was corrected a ssuming the Galactic extinction law 
of lCardelli et"aH (|l989l ), with A v = 0.085 (NGC 6868) and 
A v = 0.332 (NGC 5903) using the E(B - V) values taken 
from NED and assuming Ay = 3E(B — V). The spectra of 
NGC 6868 and NGC 5903, corrected for redshift and fore- 
ground reddening, are shown in Figs. [3] and [4] respectively. 
The distances of NGC 6868 and NGC 5903, computed as- 
suming the Hubble constant Ho = 75kms _1 , are 38Mpc 
and 34Mpc respectively. Accordingly, one arcsec is equiva- 
lent to 184 pc for NGC 6868 and 164 pc for NGC 5903. 



3 LICK INDICES 

NGC 6868 and NGC 5903 present strong absorption lines 
of neutral iron and sodium, as well as Mg2 and TiO 
bands in their visible spectra (Figs. |3] and H} that, to- 
gether with the continuum distribution, can be used to esti- 
mate the age and metallicity of the stellar po pulation. In 
the present work we adopt the Lick system (|Faber et al.l 
1 19851 ; I Worthed 1 1994h to measure equivalent widths (EWs) 
for the absorption features Mg2 A5176, FelA527o, FeL\5335, 
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Table 1. Correction factors C(<r) 
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Figure 4. Same as Fig.|3]for NGC 5903. 



FeiA5406, FeiA5709, FeiA5782, NaiA5895 and Tioa6237- Con- 
tinuum values, normalized at A5870A, were measured at 
A5300, 5546, 5650, 5800, 6173, 6620, and 6640 A. The velocity 
dispersion in NGC 6868 varies from ~ 260 ± 20 km s" 1 in the 
center to ~ 220±30kms _1 in the external regions and from 
~ 240 ± 20 kms -1 in the cent er to ~ 200 ± 40 k ms" 1 in the 
external regions of NGC 5903 l|Caon et al.ll200Gi ). Therefore, 
the observed EW values were corrected for line broadening 
due to stellar velocity dispersion. The spectrum of the G 
giant star HR5333, assumed to have zero intrinsic velocity, 
was broadened with a series of Gaussian filters with veloc- 
ity dispersion a varying from to 300 kms -1 in steps of 
10 kms" 1 . 

For each absorption feature considered we calculate an 
empirical correction index C(<r), so that C(<t)a/ 9 2 = m g2(°) ' 
and C(a)Fei = Wx(cr) • ^ e calculate the correction index 
C( a) for each extrac ted spectrum using the a values given 
bv lCaon et al.l (|2000h . The resulting correction factors, com- 
puted for each feature in all extractions, for both galax- 
ies are given in Table [1] EWs measured in the spectra of 
NGC 6868 and NGC 5903 are given in Table [3 Except for 
EW(Mg2 A5176), which is meas ured in magnitude , the re- 
maining EWs are given in A fe.g. lRickes et al.ll2004l ). Differ- 
ent values of EWs are computed allowing for three choices of 
the continuum level around a given feature, which basically 
reflect the subjectiveness associated to continuum determi- 
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nation. EWs given in Table [2]correspond to the average val- 
ues obtained from this process; the uncertainties correspond 
to the respective standard deviation. 

The radial dependence of the EWs are shown in Figs. [5] 
and© respectively for NGC 6868 and NGC 5903. Indices of 
NGC 6868 consistently present a negative gradient (Fig. [5} 
that indicates an overabundance of Fe, Mg, Na and TiO in 
the central parts with respect to the external regions. 

A conspicuous negative gradient of the indices 
Mg 2 A5176, FeiA5270, FeiA5335, NIaA5895 and TiO A 6237 is ob- 
served in NGC 5903 (Fig. |6j). 

As shown in Fig. [7] Adg2 correlates both with FeiA5270 
and FelA5335, which suggests that these elements probably 
underwent the same enrichment process in NGC 6868. How- 
ever, only a marginal correlation of Mgi and FelA5270 occurs 
in NGC 5903. 

If the galaxy gravitational potential is strong enough to 
retain the gas ejected by stars and supernovae, it eventually 
migrates to the galaxy's central regions. New generations of 
stars will be more metal rich at the center than in the exter- 
nal parts, and a steep radial negative metallicity gradient is 
established as a mass-dependent parameter. Thus, a galaxy 
formed essentially through a monolithic collapse is expected 
to present well-d efined correlations between metallicity gra - 

~ lOgando et ail 2005). 



1993a; 



die nts and mass ( Carollo et al.l 

ICarollo et al ( 1993ah investigated the correlation of the 
gradient dMg2 /dlog r with galaxy mass for a sample of 
42 elliptical galaxies, including NGC 6868 and NGC 5903. 
They found that the Mgi gradient presents a bimodal trend 
with mass, in the sense that for galaxy mass smaller than 
~ 10 11 Mq, the gradient increases with mass. No pattern in 
the gradient was observed for larger masses. The increase in 
the M<72 gradient with mass was interpreted as the result 
of a dissipative collapse associated with star formation at 
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Table 2. Equivalent widths measured in the spectra of NGC 6868 and NGC 5903 
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3.29 ±0.23 


3.21 ±0.41 


2.52 ± 0.04 


1.19 ±0.06 


1.07 ±0.02 


4.68 ±0.04 


4.52 ± 0.15 


7.32S 


0.18 


1.38 


0.23 ± 0.01 


3.57 ±0.17 


3.08 ±0.07 


2.35 ± 0.10 


1.15 ±0.04 


0.98 ±0.02 


4.54 ±0.03 


4.68 ± 0.17 


9.76S 


0.24 


1.84 
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0.46 0.26 ±0.01 3.61 ±0.13 3.01 ± 0.07 2.41 ±0.16 1.15 ±0.03 1.10 ±0.01 5.48 ± 0.01 5.77 ±0.15 



0.92 



0.22 ±0.01 3.51 ±0.16 3.06 ± 0.06 2.36 ±0.17 1.15 ±0.03 0.92 ± 0.03 4.88 ± 0.03 5.58 ±0.16 
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1.84 0.21 ±0.01 

2.30 0.21 ±0.01 
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3.22 0.17 ±0.01 



3.23 ±0.20 2.73 ±0.11 2.48 ± 0.05 0.81 ± 0.01 1.02 ± 0.08 4.59 ±0.05 5.28 ±0.19 

3.13 ±0.09 2.68 ±0.11 2.43 ±0.13 1.04 ±0.04 0.97 ±0.02 4.27 ±0.03 4.39 ±0.14 

3.01 ±0.26 2.85 ±0.04 2.21 ± 0.09 1.02 ±0.04 0.90 ± 0.07 4.17 ±0.09 4.77 ±0.23 

2.94 ±0.21 2.81 ±0.03 2.41 ±0.16 0.92 ± 0.04 1.10 ±0.13 4.10 ±0.12 4.45 ± 0.23 

2.60 ±0.15 2.36 ±0.16 1.77 ±0.13 0.77 ±0.02 0.99 ± 0.02 3.72 ±0.11 4.37 ±0.22 
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3.37 ±0.15 


2.23 ± 0.03 


1.34 ±0.02 


1.01 ±0.02 


3.76 ±0.07 


4.90 ±0.08 


24.40N 


0.66 


3.90 


0.30 ±0.02 


3.28 ±0.06 


3.33 ±0.06 


1.98 ± 0.03 


1.09 ±0.02 


0.72 ±0.10 


3.35 ±0.08 


5.48 ± 0.15 



Table Notes. Note that EW of Mg 2 A5176 is given in mag. Extraction area of all spectra of NGC6868 is 3.66(") 2 «* 0.12 kpc 2 . For 
NGC 5903 it is 4.57(") 2 ~ 0.12 kpc 2 , except for the outermost spectra, which is 9.15(") 2 ~ 0.25 kpc 2 . EWs are corrected by velocity 
dispersion and the errors are the standard deviation of three different measurements of the EW for each line. 



work in the smaller galaxies. On the other hand, the lack 
of correlation between the Mg2 gradient with mass above 
the ~ 1O 11 M0 threshold might indicate that merging of 
smaller galaxies could be the dominant formation mecha- 
nism of massive galaxies. 

To test where NGC 6868 and NGC 5903 occur in the 
dMg2/dlogr vs. galaxy mass, we compute their dynam - 
ical masses according to Ivan Dokkum fc Stanford! <|2003h . 
with the data used in the present work. The results are 
Mmgc 6868 = (3.2 ± 0.1) x 10 11 M Q and M NGC 5903 = 
(1.8 ± 0.1) x 10 11 M Q , which agree , within the u n certain - 
ties, with the values computed by ICarollo et alj (|l993al ). 
The average gradient values measured with our data 
are (dMg2/dlogr)NGC6868 = -0.07 ± 0.01 mag (") _1 and 
(dMga/dlogr) ngc 5903 = -0.03 ± 0.01 mag H" 1 . Thus 
according to the arguments developed by ICarollo et al.l 



(|l993ah . the present data are consistent with a merger origin 
for NGC 6868 and NGC 5903. 



4 STELLAR POPULATION SYNTHESIS 

The star formation history in a galaxy is a potential source 
of information not only on the age and metallicity distribu- 
tion of the stellar population components. Dating the star- 
forming episodes may provide clues to better understand 
galaxy formation and evolution processes. 

In what follows we investigate the star formation his- 
tory of NGC 6868 and NGC 590 3 by means of the stellar 
population synthesis method of iBical (|l988l ). As popula- 
tion templates we use th e synthetic star cl uster spectra o f 
iBruzual fc Chariot! (2003), built assuming a lSalpeterl (|l955T ) 
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Figure 5. Spatial distribution of EWs measured in NGC 6868. 



Figure 6. Same as Fig.[5]for the EWs of NGC 5903. 



Table 3. Age and metallicity components 



Z = 0.05 Z = 0.02 Z = 0.008 

13Gyr Al A2 A3 

5Gyr Bl B2 B3 

1 Gyr C2 



Initial Mass Function and masses in the range 0.6 ^ m s: 
120 M Q . 

Before establishing the spectral base, we note that el- 
liptical galaxies, in general, do not present recent star for- 
mation. Consequently, we only consider components older 
than 1 Gyr in the synthesis. The spectral base is made up of 
7 components with ages of 1, 5 and 13 Gyr and metallicities 
Z = 0.008, 0.02, and 0.05 (Tabled}. These metallicities are 
taken as representative of the sub-solar, solar and above- 
solar ranges. 

The stellar population synthesis provides directly the 
flux fractions (relative to the flux at A5870) that each age 
and metallicity template contributes to the observed spec- 
trum. The sum of the template spectra, according to the 
individual flux fractions, should be representative of the stel- 
lar population contribution to t he observed spectru m. A full 
description of the method is in lRickes et al.l (2004). 

As a caveat we note that the observed spectra cover 
a relatively short spectral range and, thus, contain a re- 
duced number of indices sensitive to age and metallicity 
that, in principle, could result in a small number of con- 



straints for the synthesis, especially in terms of metallicity. 
However, the sum over metallicity shows that in the cen- 
tral regions of NGC 6868 and NGC 5903, the 13 Gyr pop- 
ulation contributes with 70 ± 16% and 56 ± 12%, respec- 
tively. The 5 Gyr population contributes with 28 ± 15% and 
25 ± 12%, respectively. These results suggest the presence of 
at least two populations of different ages in both galaxies, 
especially in the central parts. The stellar population syn- 
thesis results are summarized in Fig. [8] which shows ~ la 
differences in flux contribution between the 5 and 13 Gyr 
populations, along most of the galaxy's radial extent. Inter- 
estingly, the flux contribution of the 13 Gyr population in 
NGC 6868 increases towards the central region, while that 
of the 5 Gyr decreases, within uncertainties. Flux fractions 
distribute almost symmetrically with respect to the center 
of NGC 6868. In NGC 5903, both populations present al- 
most uniform flux contributions, within uncertainties. The 
1 Gyr population presents significant flux contributions only 
in the external regions (\R\ > 2.3 kpc) of NGC 6868, and in 
the central parts of NGC 5903. 

We illustrate the stellar population synthesis in Fig. [5] 
where we superimpose on the central spectra of NGC 6868 
and NGC 5903, the respective synthesized population spec- 
tra. The residuals in the model fit to NGC 6868 blueward of 
~ 5500 A are probably due to bad calibration at the tail of 
the spectrum. The stellar population-subtracted spectrum 
of NGC 6868 present conspicuous emission lines, especially 
Hydrogen Balmer, [Nil] and [SII] . NGC 5903, on the other 
hand, has no evidence of ionized gas. 
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Figure 7. Correlations among selected Lick indices measured in 
NGC 6868 (left panels) and NGC 5903 (right panels). 



5 METALLICITY AND IONIZED GAS 
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Figure 8. Flux-fraction contribution at 5870 A of the stellar- 
population templates to the spectra of NGC 6868 (top panel) and 
NGC 5903 (bottom panel). In both cases the 5 and 13 Gyr com- 
ponents dominate the spectra, from the center to the external 
parts. 



5.1 Metallicity 

Inferences on the metallicity of NGC 6868 and NGC 5903 
can be made comparing the observed Lick indices of 
Mg 2 A5176, Fe/5270 e Fe/5335 with thos e derived from 
single-aged stellar populat ion (SSP) model s (|Thomas et al.l 
l2003l . iBuzzoni et al.lll993 ). that assume a ISalpeterl ( |l955l ) 
initial mass function and 10 Gyr of age. 

Lick indices computed from SSP models for different 
metallicities and [a/Fe] ratios together with our data are 
shown in Fig. 1101 with NGC 6868 in the top panel and 
NGC 5903 in the bottom. The central parts of NGC 6868 
(\R\ < 0.5 kpc) present a deficiency of alpha elements with 
respect to iron (—0.3 < [a/Fe] < 0.0) and an above- 
solar metallicity ([Z/Zq] « +0.3. The external parts, on 
the other hand, present a higher, roughly uniform distribu- 
tion of ratios ([a/Fe] ps +0.3) and sub-solar metallicities 
([Z/Zq] « —0.33). The overabundance of the [a/Fe] ratios 
has been interpreted as a consequence of the chemical en- 
richment pro duced by type II supernovae with respect to 
type la ones (jldiart et al.ll2003l ). Thus, differences measured 
in [a/Fe] imply different star- formation histories for the cen- 
tral and external regions of NGC 6868, which is consistent 
with the results derived from the stellar population synthesis 
(Sect. 2]). A similar conclusion applies to NGC 5903 (bottom 
panel), since the central parts have metallicities in excess of 
[Z/Zq] ss +0.35 and [a/Fe] fa 0.0, while the external parts 



have metallicities between solar and [Z/Zq] w +0.35 and 
present an enhanced ratio [a/Fe] ~ +0.3. 



5.2 Ionized gas 

Emission gas has be en recently detected in a large number 
of elliptical galaxies (jPhillips et alJll986h . However, its ori- 
gin and the nature of the ionization source have not yet 
been conclusively established. NGC 6868 presents conspic- 
uous emission lines not only in the central region, but in 
spectra extracted up to R ~ 17" ~ 3.1 kpc (Fig. In 
what follows we investigate properties of the emission gas 
using fluxes of the lines Ha, [ATI] aa 6548,6584, [O/]a6300 
and [S//]aa 6717,6731 > measured in the stellar population-free 
spectra, i.e., those resulting from the subtraction of the re- 
spective population templates (Sect. |4}. The emission line 
fluxes were measured fitting Gaussian to the profiles. 

Flux ratios with respect to Ha are given in Tab. [4] and 
the spatial distribution of the measured ratios tfjJ, ^jtzt 
and are shown in the left panels of Fig. 1111 

Assuming case B recombination, the number of ioniz- 
ing photons can be computed from H a luminosity using 
the equation Q(H) = igf "^(h^t) <|Osterbrockl Il989l ). 
where as (H°,T) is the total recombination coefficient, and 
cibo (H ,T) is the recombination coefficient for Ha. Values 
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Figure 9. Stellar population synthesis of NGC 6868 (top panel) 
and NGC 5903 (bottom panel). Notice that the subtraction of the 
stellar population enhanced Ha in emission in NGC f 



Table 4. Emission line parameters of NGC 6868 



R{") 
(1) 


Ha 
(2) 


[Nil] 
Ha 

(3) 


[ST/] 
Ha 

(4) 


ion 

Ha 

(5) 


0.00 


3.36 ±0.08 


2.81 ±0.21 


1.23 ±0.10 


0.51 ±0.04 


2.44S 


3.39 ±0.10 


2.90 ±0.91 


1.18 ±0.12 


0.33 ±0.04 


4.88S 


1.81 ±0.11 


2.46 ±0.21 


0.92 ±0.11 


0.51 ±0.07 


7.32S 


1.20 ±0.18 


2.36 ±0.36 


0.88 ±0.40 


0.92 ±0.17 


9.76S 


0.82 ±0.21 


2.36 ±0.61 


1.20 ±0.31 


1.23 ±0.34 


2.44N 


3.83 ±0.13 


2.92 ±0.32 


1.32 ±0.11 


0.29 ±0.03 


4.88N 


2.86 ±0.10 


2.65 ±0.22 


1.30 ±0.12 


0.26 ±0.04 


7.32N 


2.34 ±0.13 


2.87 ±0.64 


1.25 ±0.10 


0.19 ±0.05 


9.76N 


2.08 ±0.16 


3.11 ±0.42 


1.39 ±0.11 


0.35 ± 0.06 


12.20N 


1.49 ±0.13 


2.88 ±0.33 


1.47 ±0.41 


0.56 ±0.11 


14.64N 


1.39 ±0.18 


2.56 ±0.57 


1.23 ±0.17 




17.08N 


1.06 ±0.20 


2.43 ±0.48 


1.20 ±0.23 





Table Notes. Ha flux in col. 2 is given in 10 ergs cm . 

of Q(H) and L H a for NGC 6868 are given in Tab. El and 
their spatial distributions are shown in panels (d) and (e) of 
Fig. ED 

Fig. [TT] also shows that in all extracted spectra of 
NGC 6868, the ratio [iY - fJ '* 6584 is larger than 1 and [S "' A6731 

Ha ° Ha 

is smaller than 1, within uncertainties. 

With the above facts in mind we consider three possi- 
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Figure 10. Selected Lick indices measured in NGC 6868 (top 
panel) and NGC 5903 (bottom panel) are compared to those com- 
puted from SSP models for a range of metallicity and [a/Fc] ra- 
tios. Ellipses and boxes indicate the central and external regions 
respectively. 



Table 5. Ha luminosity and number of ionizing photons in 
NGC 6868 



R 

(") 



(10 38 ergs^ 1 ) 



Q(H) 
(lO^s- 1 ) 



0.00 


5.8 ± 1.2 


2.6 ±0.5 


2.44S 


5.9 ± 1.2 


2.6 ±0.5 


4.88S 


3.1 ± 0.6 


1.4 ±0.3 


7.32S 


2.1 ± 0.4 


0.9 ±0.2 


9.76S 


1.4 ±0.3 


0.6 ±0.1 


2.44N 


6.6 ± 1.3 


2.9 ±0.6 


4.88N 


4.9 ± 1.0 


2.2 ±0.4 


7.32N 


4.0 ±0.8 


1.8 ±0.4 


9.76N 


3.6 ± 0.7 


1.6 ±0.3 


12.20N 


2.6 ± 0.5 


1.1±0.2 


14.64N 


2.4 ±0.3 


1.1 ±0.1 


17.08N 


1.8 ± 0.4 


0.8 ±0.2 



Table Notes. Ha luminosity and number of ionizing photons 
computed for the extractions given in col. 1 of NGC 6868. 



ble scenarios to infer the nature of the ionization source in 
NGC 6868. They involve the presence of star clusters, post- 
AGB stars and an AGN. 
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Table 6. Hypothetical ionizing star cluster 
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Figure 11. Spatial variation of the ratios 



[oih 



^ 30 ° (panel a), 

jjr* 6584 (b) and |jj ^ a q 6731 (c). Panels (d) and (e) show the 
spatial distribution of the number of ionization photons and Ha 
luminosity, respectively. The dashed line in panels (a) and (c) 
discriminates between different ionizing sources (Sect. IB",2.2I I. 



5.2.1 Star cluster 

To test this hypothesis we build a star cluster that pro- 
duces a total num ber of io n izing photons of the order of 
10 51 s" 1 following ISalpeterl (] 19551 ) initial mass function, 
4>{m) oc m~' 1+x ', with \ ~ 1-35, considering stars in the 
mass range 0.1 — 100 Mq. For the relation of number of 
ionizing photons (Nl v ) with st ellar mass we use the li- 
brary of stellar atmospheres of lKuruc3 (1 19791 ). In cols. 2 
and 3 of Table [5] we provide Nl v for ionizing stars in the 
mass range 10 — 30 Mq . In this mass range, the number 
of ionizing photons can be related to stellar mass by the 

~i.9 



39 ( m V 



Thus, 



approximation NL y (m) 3.9 x 10 
the total number of ionizing photons can be computed from 
Jio N Ly {m)<f){m)dm. 

The resulting distribution of stars in terms of spectral 
type of the hypothetical ionizing star cluster is given in col. 4 
of Tab. [6] for representative spectral types. Because of the 
large number of 05 stars, such a star cluster would be very 
young, luminous and have a mass of ~ 10 6 Mq . 

We conclude that ionization by such a massive star clus- 
ter can be ruled out, because with more than 200 05 stars 
it should be bright enough to be easily identified in images. 
Besides, its flux contribution to the observed spectra should 



Spectral type 


rn 




Number of stars 




(Mq) 


(10 48 s" 1 ) 




(1) 


(2) 


(3) 


(4) 


05 


30 


2.00 


226 




22 


0.25 


470 


B0 


17 


0.082 


860 




15 


0.043 


1154 




13 


0.015 


1616 




10 


0.0021 


2993 



be conspicuously detected by the stellar population synthe- 
sis (Sect. g|| . 



5.2.2 Post-AGB stars 



Since a typical Post-AGB star emits ~ 10 47 s _1 (|Binettd 
Il994l ). about 10 4 such stars would be necessary to produce 
the amount of ionizing photons measured in each extracted 
region (Tab. [fa]). All spectra of NGC 6868 were extracted 
from regions with an area of ~ 0.19 kpc 2 . This implies a pro- 
jected number density of Post-AGB stars of ~ 8 x 10 4 kpc -2 
scattered throughout th e R » 3.3 kpc central region. Such 
a density is not absurd l|Binettell 19941 ). however, Post-AGB 
star ionization models for a variety of ionization parame- 
ters predict that the intensity ratio [O/]a6300 /Ha should 
be smaller than unity, while those of [S I T]a673i/ Ha and 
[A r //]A6584/ 'Ha should be smaller than w 1.5. The first two 
conditions are met by the measurements of NGC 6868, but 
the ratios [NII]\e5si/Ha are all above 1.5 (Fig |ll[l . 

The presence of a significant number of post-AGB stars 
in early-type galaxies p r oduce s a conspicuous excess in the 
UV flux, e.g. lBica et atl (ll996l). NGC 6 8 68 is included in the 
galaxies investigated by iBonatto et all (|l996l ). among those 
with a flat UV spectrum, which is consistent with the above 
arguments against a numerous population of post-AGB stars 
in this galaxy. 

5.2.3 Active Nucleus 

To test this hypothesis we build diagnostic-diagrams 
involving the ratios [NII]\e5S4/Ha, [07]a630o/ Ha and 
[5//]A67i7,3i/^ a , using as comparison the values measured 
in a sample if nearby galaxies wi th low luminosi ty (Lua < 
2 x 10 39 ergs _1 ) active nuc l eus jHo et al.|[l997h . Here we 
use the galaxies in iHo et all (1 19971 ) sample with unambigu- 
ous classification as LINER in spiral, LINER in elliptical, 
Seyfert or Starburst. The restricted sample amounts to 128 
galaxies. 

The results are shown in Fig. 1121 where in both panels 
we see that the values measured in NGC 6868 consistently 
fall in the locus occupied predominantly by LINERs. 

These results, together with those of the stellar popu- 
lation synthesis, suggest that the main source of gas ioniza- 
tion in NGC 6868 is non-thermal. However, as a caveat we 
note that shocks can also produce emission-line ratios sim- 
ilar to those measured in NGC 6868. In this sense, the ra- 
tios [NII]xe5S4/Ha and [SII]xsnr, 31 /Ha alone cannot dis- 
criminate betw een the AGN and shock s as t he ionization 
source. Indeed. [storchi-Bergmann et al.l (|l996l ) show that a 
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Figure 12. Emiss ion-line diagnost ic diagrams with a subsample 
of the galaxies in iHo et aD l|l997t ) plotted as comparison. The 
locus of NGC 6868 in both diagrams is consistent with that of 
LINERs. 



with 70 ± 16% and 56 ± 12%, respectively. The 5 Gyr popu- 
lation contributes with 28 ± 15% and 25 ± 12%, respectively 
in flux fraction at 5870 A. 

The central regions of NGC 6868 and NGC 5903 present 
a deficiency of alpha elements with respect to the exter- 
nal parts, which suggests different star-formation histories 
in both regions. With respect to the metallicity, the cen- 
tral regions of both galaxies have higher [Z/Zq] values than 
the external parts. The range in metallicity spanned by the 
sampled regions of NGC 6868 appears to be larger than in 
NGC 5903. 

Concerning the ionized gas in NGC 6868, the ratios 
[A''77]A6584/-H'a, [OI]\630o/Ha and [SI 7]a6717,3i/ Ha consis- 
tently suggest the presence of a LINER at the galaxy center. 
These results, together with the stellar population synthesis, 
suggest that the main source of gas ionization in NGC 6868 
is non-thermal, produced by a low-luminosity active nucleus, 
probably with some contribution of shocks to explain ion- 
ization at distances of ~ 3.5 kpc from the nucleus. 
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low-luminosity AGN cannot ionize gas up to distances of 
~ 3.5 kpc. Besides, they suggest that ionization in the ex- 
ternal region is produced by shocks, which accounts for the 
constancy of the [A r 77]A6584/^ a and [S , 77]a6717,3i/7-7o: ra- 
tios. 



6 DISCUSSION AND CONCLUDING 
REMARKS 

The stellar population, metallicity distribution and ionized 
gas in NGC 6868 and NGC 5903 have been investigated in 
this paper by means of long-slit spectroscopy and stellar 
population synthesis. Lick indices of both galaxies consis- 
tently present a negative gradient that indicates an over- 
abundance of Fe, Mg, Na and TiO in the central parts 
with respect to the external regions. We found that Mg2 
correlates both with Fel,\5270 and Fel,\5335, which suggests 
that these elements probably underwent the same enrich- 
ment process in NGC 6868. However, in NGC 5903 only a 
marginal correlation of Mg2 and FelA5270 occurs. 

Galaxy mass and the Mg2 gradient computed in the 
present work are consistent with previous results that sug- 
gest that NGC 6868 and NGC 5903 were formed by merger 
events. In addition, the stellar population synthesis clearly 
shows the presence of at least two populations of different 
ages in both galaxies. Particularly in the central regions of 
NGC 6868 and NGC 5903, the 13 Gyr population contributes 
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